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Abstract 



Using connection with quantum field theory, the infinitesimal covariant abelian gauge 
transformation laws of relativistic two-particle constraint theory wave functions and po- 
tentials are established and weak invariance of the corresponding wave equations shown. 
Because of the three-dimensional projection operation, these transformation laws are in- 
teraction dependent. Simplifications occur for local potentials, which result, in each formal 
order of perturbation theory, from the infra-red leading effects of multiphoton exchange 
diagrams. In this case, the finite gauge transformation can explicitly be represented, with 
a suitable approximation and up to a multiplicative factor, by a momentum dependent 
unitary operator that acts in x-space as a local dilatation operator. The latter is utilized 
to reconstruct from the Feynman gauge the potentials in other linear covariant gauges. 
The resulting effective potential of the final Pauli-Schrodinger type eigenvalue equation 
has the gauge invariant attractive singularity a^/r^, leading to a gauge invariant critical 
coupling constant etc = 1/2. 

PACS numbers: 03.65.Pm, ll.lO.St, 12.20.Ds, ll.80.Fv. 

Keywords: Relativistic wave equations, Bethe-Salpeter equation, constraint theory, 
gauge transformations. 
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1 Introduction 



The knowledge of the behavior of Green's functions under gauge transformations in QED 

0, plays a crucial role in proving gauge invariance of observables and in particular 
of bound state energies @, H- In practical calculations, however, one generally uses 
approximations to exact equations, which necessitate a close control of the degree of 
realization of the various symmetries of the system under study. 

In QED, it has appeared that the Coulomb gauge is the most convenient gauge for 
treating the bound state problem, since it allows the optimal expansion of the Bethe- 
Salpeter equation P, |^, || around the nonrelativistic theory 0, |], |10|, |Tl], [12|. Covariant 
gauges produce, at each formal order of perturbation theory, spurious infra-red singu- 
larities that are cancelled only by higher order diagrams |TB|, |T^ and thus become of 
less practical interest. The main disadvantage of the Coulomb gauge is, however, its 
noncovariant nature, which does not allow its incorporation in covariant equations. 

From this viewpoint, constraint theory, which leads to a manifestly covariant three- 
dimensional description of two-body systems |16], [l^, has opened a new perspective in 
the subject. It was shown |T^ , that the expansion of the Bethe-Salpeter equation 

around the constraint theory wave equations in the Feynman gauge (as well as for scalar 
interactions) is free of the abovementioned diseases of covariant gauges and allows a 
systematic study of infra-red leading effects of multiphoton exchange diagrams; the latter 
can then be represented in three-dimensional x-space as local potentials. Summing the 
series of these leading terms one obtains a local potential in compact form ||20|, which 
is well suited for a continuation to the strong coupling domain of the theory or for a 
generalization to other effective interactions. 

The purpose of the present paper is to investigate the forms of the local potentials 
in linear covariant gauges and to find the relationships between them. Our approach 
is accomplished in two steps. First, using the connection of the constraint theory wave 
equations with the Bethe-Salpeter equation, we determine the infinitesimal gauge trans- 
formation laws of constraint theory wave functions and potentials for fermion-antifermion 
systems and establish the (weak) invariance of the corresponding wave equations. Second, 
specializing to the local potential approximation, we show that the finite gauge transfor- 
mation can be represented, with a suitable approximation and up to a multiplicative 
factor, by a momentum dependent unitary operator that acts in x-space as a local di- 
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latation operator. Under the action of the latter, the potentials undergo two kinds of 
modification. In the first one, they are changed in a form invariant way, by the replace- 
ment of their argument r (the cm. interparticle distance) by a function r(^), where ^ 
is the gauge variation parameter. Apart from a rapid variation near the origin, r(^) is 
essentially dominated by its large-distance behavior, in which r is simply shifted by a con- 
stant value. In the second kind of modification, certain parts of the spacelike components 
of the electromagnetic potential are functionally modified. The combination of these two 
types of modification allows, in particular, the reconstruction of the potentials in linear 
covariant gauges from the knowledge of the potential in the Feynman gauge. 

Studying the short-distance behavior of these (three-dimensional) potentials we find 
that the effective potential that appears in the final Pauli-Schrodinger type eigenvalue 
equation has the gauge-invariant attractive singularity o^jr^. As is known, to such a 
singularity there corresponds a critical value etc of the coupling constant, — 1/2, at 
which value the fall to the center phenomenon occurs. 

Constraint theory thus provides us with a three-dimensional framework in which the 
local potential approximation consistently fulfils the requirement of gauge invariance of 
the theory. These results, while deduced in the framework of QED, might also survive, 
with appropriate adaptations, to the incorporation or consideration of other types of 
interaction. 

The plan of the paper is the following. In Sec. 2, the gauge transformation properties 
of the Green's functions and of the Bethe-Salpeter wave function are reviewed. In Sec. 3, 
the infinitesimal gauge transformation laws of the constraint theory wave functions and 
potentials are determined. The case of the local approximation of potentials is considered 
in Sec. 4. The general properties of the gauge transformations in the local approximation 
are studied in Sec. 5. A summary and concluding remarks follow in Sec. 6. 
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2 Gauge transformations of Green's functions and 
wave functions 



Gauge transformation laws of Green's functions in QED can be obtained with several 
equivalent methods: i) by considering the operator changes in the field operators ii) 
by modifying in the functional integral the gauge fixing condition 0]; iii) by using the 
Ward-Takahashi identities . In the present work we shall consider only linear covariant 
gauges, characterized by a parameter ^ ; the photon propagator is then: 

= (2.1) 

The reference gauge is taken here to be the Feynman gauge and the transformations 
express a Green's function calculated in the gauge ^ with respect to its value in the 
Feynman gauge = 0). More generally, these transformations concern the passage from 
a gauge to a gauge ^2, with ^ = ^2 - 6- 

The transformation law of the unrenormalized Green's function of a charged particle 
(a boson or a fermion) is: 

G^{x) = exp{ie^^{A{x) ~ A{0))} G{x) , (2.2) 

where e is the unrenormalized electric charge of the particle and A is defined as: 

Notice that the transformation law is spin independendent. 

For a two-particle Green's function, with particle 1 representing a fermion (boson) 
with charge ei and particle 2 an antifermion (boson) with charge —62, the transformation 
law is [|, ig: 

G^ixi, X2; 2/1, 2/2) = exp{ze[e?(A(xi - yi) - A(0)) + el{A{x2 - 2/2)) - A(0))]} 

X exp{z^eie2[A(a;i - X2) + A(?/i - 2/2) - A(a;i - 2/2) - A(yi - X2)]} 
xG{xi,X2;yi,y2) ■ (2.4) 



One also deduces from Eq. (|2.4| ) the gauge transformation law of the Bethe-Salpeter 
wave function. We assume that the two-particle system is neutral in charge, hence: 

ei = 62 = e . (2.5) 



The time separation between outgoing and ingoing states is defined as T = ((x? + X2) — 
iUi + vV))/'^- -^o^ large time separations (T 00) and finite values of x\ — x\ and y1 — 1/2, 
G satisfies the cluster decomposition |^: 

G(xi,X2;?/i,?/2) = vJ2^n(^i^^2)^niyi,y2) , (2.6) 

where t] = +1 for bosonic fields and 77 = —1 for fermionic fields; the sum over n cor- 
responds to a complete set of states; $„ is a generalized Bethe-Salpeter wave function 
with respect to the intermediate state \n > and becomes the usual Bethe-Salpeter wave 
function when |r;, > is a bound state P, 0. Although the function A has a logarithmic 
increase at infinity, the charge neutrality condition (|2.5|) leads to a mutual cancellation of 
singular terms and one obtains |r5| : 

X2) = e^e'^(^(^i - ^2) - A(0)) X2) . (2.7) 

The infinite factor A(0) in Eqs. ( p.2|) , (|2.4] ) and ( p.7|) determines the transformation 
law of the wave function renormalization constant ^'2 [0. In the remaining part of this 
work, we shall not consider radiative corrections and hence shall ignore gauge transfor- 
mation effects coming from photons associated entirely with one particle; only exchanged 
photon contributions will be taken into account. [More generally, one may consider dif- 
ferent gauges for exchanged photons and for photons entering in radiative corrections 
0.] Except in particular limiting procedures, such as in Eq. (|2.6| ), there are no, in gen- 
eral, cancellations between the two kinds of contributions, since they are concerned with 
different variables. 

While the behavior under gauge transformations of Green's functions and wave func- 
tions is very simple (in x-space), the same is not true for the inverses of Green's functions, 
the vertex functions and the scattering amplitudes. In this case, only infinitesimal gauge 
transformations have relatively simple forms. They can be obtained either by inverting 
the Green's functions or by using the Ward-Takahashi identities. 

The infinitesimal change of the inverse of the four-point Green's function is: 

6^G-^ = -G-\S^G)G-^ . (2.8) 



Equations (|2.8|) and (|2.4|) can then be used to show the invariance, in a weak form, of the 



Bethe-Salpeter equation. The latter can directly be written in terms of G ^ as: 

G-^<!> = . (2.9) 



Invariance of this equation under gauge transformations requires: 

{6^G-^)<I> + G-^6^<l> = G-U-{6^G)G-^<I> + 6^<^) = 0. (2.10) 



Replacing in this equation S^G and 5^$ by their expressions ( |2.4| ) and (p.7|), respectively, 
one obtains: 

A{xi - X2) j d'^zid'^Z2G~^{xi, X2; zi, Z2)^{zi, Z2) 

— J d'^x[d'^x'2d'^yid'^y2d'^Zid'^Z2G~^{xi, X2; x[, x'2){/S.{x[ — ^2) + ^(2^2 ~ Vi)) 
xG{x[,X2;yi,y2)G-\yi,y2;zi,Z2) ^{zi,Z2) = 0. (2.11) 

(The contribution of 5^$ has been cancelled by one of the terms of 6^G.) The first term 
vanishes on account of Eq. ( p.9|) . In the second and third terms, the gauge propagators A 
acting multiplicatively on G, join one of the x"s to the y of the other line; in momentum 
space, they produce crossed diagrams with G; these do not have bound state poles at 
the positions of the direct diagram poles; therefore, the two products AG cannot prevent 
the factor from vanishing, on account of Eq. (|2.9| ). One thus establishes the weak 

gauge invariance of the Bethe-Salpeter equation. 

This conclusion is not changed when one includes radiative corrections relative to 
fermion or boson lines and vertices. In particular, radiative corrections that join x' to y 
on the same line, imply in momentum space integrations on the bound state pole position, 
which is then transformed into a cut. 

To display the gauge transformation property of the off-mass shell scattering ampli- 
tude, we pass into momentum space. The infinitesimal form of Eq. ( p.4|) is: 



f d'^k r 

5^G(pi,-p2]p'i,-p'2) = ie^^i j j^--^A{k) G{pi - k, -{p2 + k);p[, -P2) 
+G{pi, -P2,p'i + k, -(p2 - A;)) - G{pi - k, -p2]p'i, -{p'2 - k)) 
-G{pi, -{p2 + k);p[ + k, -P2) 



(2.12) 



This equation is graphically represented in Fig. 1. We notice that if G has a (simple) 
bound state pole, then in the right-hand side of Eq. ( |2.12| ) the first two terms have also the 
same pole. However, this pole being simple, one concludes that the gauge transformation 
does not change the bound state energy |]. (Otherwise, the singularity structure in 
the right-hand side of Eq. ( |2.12| ) would be that of a double pole.) 
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The scattering amplitude T is defined from tlie Green's function G as: 



T — ^ [G* — GqJG-i^ ^G 



2 5 



(2.13) 



where Gi and G2 are the external particle propagators and Go is their product: 



Gq — G1G2 ■ 



(2.14) 



The infinitesimal gauge transformation law of T is then: 



G-,\p,) - G-,\p\) 



G,\ 



-P2) 



G,H 



-P2) 



(2^ 



G^\pi)G,{p, - k)T{p, - k, -{p2 + k);p[, -p'2) 



ip^ + k))G-^'' 



^2\~yP2 -r l^))'^2 \-P2) 

+G^\-p'^)G,{-{p', - k))T{p,, -p,;p[ + k, -ip', - k))GM + k)G^\p[) 
-Gi\p,)G,\-p',)G,{p, - k)G2{-{p', - k))T{p, - k, -p2;p[, -{p', - k)) 



-TiPi, -{P2 + ky,p[ + k, -p',)GM + k)G2{-{p2 + k))G^\p[)G^^ 



-P2) 



(2.15) 



Taking the external particles on their mass-shell, one immediately deduces from this 
equation gauge invariance of the on-mass shell scattering amplitude. 

The unitarity of the gauge transformation can be checked with the invariance of the 
norm of the Bethe-Salpeter wave function. We first notice that, according to Eq. ( |2.4| ), 
the adjoint $ of the Bethe-Salpeter wave function [Eq. ( p.6|) ] transforms in the same way 
as $ [Eq. ( p.7|) ] (with the same sign in the exponential function); this is due to the fact 



that the gauge function A [Eq. ( p. 3D ] is imaginary for euclidean variables in x-space. The 
norm of the internal part of the Bethe-Salpeter wave function is [Q : 

dG~^ , 



ds 



—ITj 



(2.16) 



where s = P"^ = {pi + P2Y and rj was defined after Eq. (|2 

Using infinitesimal gauge transformations together with Eqs. ( |2.7| ), ( p.8| ) and ( p.l2| ) 
and restricting ourselves to hermitian kernels and time-reversal invariant interactions, we 
find (in compact notation): 
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-. dG- 



-T 9G-^ . . 
+ / A0 



ds 



ds J ds 
(AG) + (GA) - (AG) 



(GA), 



0G 



-1 



(AG) + (GA) - (AG)„. - (GA)„, 



OS 



(AG) + (GA) - (AG),, - (GA), 



ds 
G-\ 



(2.17) 



where the subscript "cr" designates the terms corresponding to the crossed diagrams [Eq. 
(|2.12|) and Fig. 1]. In the above equation, one is entitled to use the wave equation (|2^ ), 
or its adjoint, (f)G~^ = 0, as long as G^^ is not multiplied by terms having a pole at the 
bound state position in the s-channel. Therefore, the crossed diagram terms disappear 
and after some algebraic cancellations one remains with the following terms, which also 
ultimately vanish: 



+ 



~ (GA)G- V + / 0G-i(AG) 



ds 

-ze'^ej / </)G-^^(^AGG-i)0 + j 0^ (^G-^GA^G^V 



ds 



AG + GA 



. 



(2.18) 



Let us finally remark that the gauge propagator A, having a singular infra-red be- 
havior, may lead, during the evaluation of Feynman diagram integrals, to infra-red diver- 



gences. However, the integrals involved in the variations 6^G [Eq. ( 2.12 )] and 6(T [Eq. 
( |2.15| )] are actually globally infra-red regular, although individual terms are divergent, as 
can be checked by taking in the integrands the limit A; ^ 0. In later calculations, some of 
the terms present in 5^G or 5(1 will be dropped, for they do not contribute to the bound 
state poles. It should be understood, in such cases, that their infra-red divergent parts 
are kept, in order to maintain the regularity of the remaining integral. 
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3 Gauge transformations of constraint theory wave 
functions and potentials 



Constraint theory [T^, |T6|, 0] allows the elimination of the relative energy variable of the 



two particles by means of a manifestly covariant equation. The choice of the latter is not 
unique and generally various choices are related one to the other by canonical or wave 
function transformations. We choose the following constraint: 

[ (pI - pI) - {ml - ml) ] ^(xi, X2) = , (3.1) 

where pi and p2 are the momentum operators of particles 1 and 2, mi and m2 their respec- 
tive masses and \E' is the constraint theory wave function. We use standard definitions for 
the total and relative variables: P = pi +P2, p = {pi ~P2)/2, x = Xi —X2, M = mi +m2. 
For states that are eigenstates of the total momentum P, we define transverse and lon- 
gitudinal components of four-vectors with respect to P and denote them with indices T 
and L, respectively; thus: xT^ = — x.PP^, xl = x.P, Pfj, = P^j^fP^, Pl = ^/P^, 



r = ^ —x^"^. 

With these definitions, constraint ( |3.1| ) can be written in the following form: 

C(p) = 2PlPl - {ml - m^) ^ . (3.2) 

Constraint C allows the elimination of the relative longitudinal momentum p^ and fixes 
the evolution law with respect to its conjugate variable xl. The internal dynamics of the 
system becomes three-dimensional, apart from the spin degrees of freedom, expressed in 
terms of the transverse vector x'^ . 

In the presence of constraint C, the individual Klein-Gordon operators become equal: 

+ + (3.3) 



^0 = {Pi - "^1 



= {P2-m2) 

C 



c 4 2^ ' 4P2 
The wave equations for two spin-0 particle systems are: 

{pl-ml-V)^!{xi,X2) = Q , (3.4a) 
{pl-ml-V)^{xi,X2) = Q , (3.4b) 

while for spin-i fermion-antifermion systems they are |]T^ pT| : 



(71 .pi - mi)^ = (-72.P2 + m2)V-^ , (3.5a) 
-72-P2 - m2)^ = (71 .pi + mi)V^ . (3.5b) 
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[The antifermion Dirac matrices 72 act on \E', represented as a 4 x 4 matrix, from the right, 
in the reverse order of their appearance.] Potentials V are Poincare invariant operators 
and depend on x only through the transverse vector x^. 

Equations ( p.4a]) -( p.4b|) and (|3.5a|) - (p.5b| ) can also be written in a unified form. Intro- 



ducing the individual particle propagators Gi and G2 (with i-factors in the numerators 
in momentum space) and their product Gq [Eq. ( ^.14| )] and defining 



go = HqG, 

Eqs. (|3^)-(|3^ and (p3a|)-(|33B[) take the form: 



(3.6) 

C{P) 



{g^' + V)^ = 0. (3.7) 

In order to establish the connection of constraint theory wave equations and potentials 
with the Bethe-Salpeter equation and related quantities, one projects, with an appropriate 
weight factor, which is chosen here to be Hq [Eq. (|3.3| )], Green's functions, scattering 
amplitudes and wave functions on the constraint hypersurface (^). Thus, defining the 
left-projected Green's function G, 

G{P,p,p') = -2zn6{Gip))HoGiP,p,p') , (3.8) 

one can iterate, in the right-hand side, G around G, repeatedly using its integral equation. 
One ends up with an integral equation satisfied by G, the kernel of which is related to 
the Bethe-Salpeter kernel K Defining the constraint theory wave function from 



the residue of G at a bound state pole, 

= 2n2PL6{G)^ = -2i7r6{G)Ho^ , (3.9) 
one finds that the wave function satisfies Eq. ( |3.7| ) with V related to the scattering 



amplitude by means of a Lippmann-Schwinger-quasipotential type equation [^, [2^] : 

V = f{l-g,f)-' , (3.10) 

(3.11) 



T(p,p') = i^np^p'] 



C(p),C{p') 

[In T, the total four-momentum conservation factor (27r) ^5'^ (P — P') has been amputated.] 
Conversely, the Bethe-Salpeter wave function $ can be reconstructed from \E'c' with 
the equation 

$ = GoT{l-gTy'^c , 9 = 2t7r6{C)go , (3.12) 
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Go and go being defined in Eqs. ( ^.141 ) and 

Equation ( p.lO| ) is the basis for the calculation of the potential of constraint theory 
from the scattering amplitude. Iterating Eq. ( |3.1CI|) with respect to T, one finds that V 
receives contributions, in addition to those of the ordinary Feynman diagrams, from new 
diagrams having at least one three-dimensional box sub-diagram, corresponding to the 
presence of the constraint factor g^. These diagrams, which we call "constraint diagrams" , 



play a crucial role in the cancellation mechanism of spurious infra-red singularities |T9|, |20 
and in the reorganiszation of the perturbation series. 

The infinitesimal gauge transformation law of the constraint theory wave function is 
obtained by starting from Eqs. (|3^ ) and using Eqs. ( p. 121) and (|2.7| ). One has: 



2P, 



c 



2Pj 



2P, 



c 



(3.13) 



where the integration inside the term {AGqT) is four- dimensional. 

The transformation law of the potential is obtained by starting from the relationship 
between the wave equation operator {g^^ + V) [Eq. ( p.7| )] and the scattering amplitude 
f [Eqs. (iAg)-(|Tl]: 



9o 



[l-Tgo 



One finds the relation 



sd9o' + v) = + v)UkT)u%' + y) 



(3.14) 



(3.15) 



which, according to the transformation law (|2.15| ), becomes 



+ (A + TGoA)Go ' - (crossed)|^o(^o"' + V) ■ (3-16) 

When this equation is applied on the wave function the operator {gQ^ + V) of 
the utmost right gives zero [Eq. ( p.TI) ] provided it is not multiplied by terms having the 
bound state pole at the same position. Therfore, the crossed terms, as well as the single 
A terms, disappear from the equation. Furthermore, in the second term of the right-hand 
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side of the equation the product {i/{2PL)){%^ + V)gQT is , according to Eqs. ( ^.141 ), the 
four-dimensional continuation (in relative longitudinal momentum) of (1 — TgQ)'~^T = V, 
which does not have a pole; hence, it does not contribute when Eq. (|3.16|) is applied on 
^f. One thus obtains the equation: 

M9o' + V)]^ = 7^^e'6a%' + V^)^^o(AGor)(l - gofy' ^ , (3.17) 

where Eq. ( ^.61 ) was used. 

The combination of the two transformation laws ( |3.13| ) and ( |3.17| ) leads to the weak 
invariance of the wave equation (|3.7|): 



6di9o' + V)<il] ^ 0. (3.18) 

The invariance of the norm of \^ can be shown in a similar way as for the Bethe-Salpeter 
wave function [Eqs. (|2l6| ) - (^) ]. The norm of the internal part of \E' is |jlO|, |l^, 

(^,V^) = -vj d^x'^AP' f ^Jg,' + V]i^ = 2Pl, (3.19) 

where is the adjoint of i/j, equal to t/j* in the bosonic case and to ['yiLl2Li^]^ in the 
fermionic case; s = P^ and rj was defined after Eq. (|2.6|) . Using the transformation laws 
( |3.13| ) and ( ^.16 ) and arguments similar to those used in Eqs. (2.1tj)-( 2.18| ) one shows the 



(weak) invariance of the norm ( p.l9| ): 

5S,i^) = 0. (3.20) 



Contrary to the four-dimensional case, the gauge transformation laws ( 3.13 ) and ( 3.17 ) 
of the three-dimensional theory explicitly depend on the interaction, a feature that ren- 
ders their evaluations rather tricky. Furthermore, the amplitude T that appears in the 
transformation is not submitted on its left to the constraint ( |3.2| ). Therefore, the quan- 
tity {i/{2PL))T{l — goT)~^ represents a four-dimensional continuation of potential V [Eq. 



(|3.1CI|) ]. Its evaluation, in the general case, cannot be done in compact form. However, 



when the local approximation is used for the potentials, simplifications occur. This case 
is considered in Sec. 4. 



13 



4 Gauge transformations in the local 
approximation of potentials 



In lowest order of perturbation theory (one photon exchange), relationship (|3.10| ) provides 
a local expression for V in three-dimensional x-space (with respect to x'^). It turns out 
that this property can also be maintained, in a certain approximation, in higher orders. 



The perturbation theory calculations effected in Ref. |20| have shown that, in the Feynman 



gauge, the infra-red leading part of the contribution of the n-photon exchange diagrams 



can be represented in (three-dimensional) a;-space by a local function of r (= y/—x^), of 
the type (a/r)^, a being the fine structure constant. The sum of these leading terms also 
yields a local function for V. Therefore, the local approximation of V can be considered 
as a sensible one: it includes not only lowest order effects, but also leading effects of 
multiphoton exchange diagrams. The use of this approximation considerably simplifies 
the resolution of wave equations, where now standard methods of quantum mechanics can 
be applied. In the rest of this article we shall limit ourselves to this approximation and 
shall correspondingly consider the transformation laws obtained in Sec. 3. 

In the fermionic case, to ensure positivity of the norm, potential V must satisfy an 



inequality ^T|]; for local potentials that commute with 7il72L; the parametrization |^| 

V = tanhl^ , (4.1) 
satisfies this condition. Introducing the wave function transformation 

^ = (coshl^)^ , (4.2) 
the norm of the internal part of the new wave function \E' becomes (in the cm. frame): 

dV 



1 + 47io72o-Po^ 



Qp2 



^} = 2Po, (4.3) 



which, for cm. energy independent potentials, reduces to the conventional free norm of 
states. 

Considering local potentials V (in x"^) composed of combinations of scalar, pseu- 
doscalar and vector potentials, 

T T 

V = V, + 715725^3 + l^l2i9iuV2 + gllU, + ^T4) , (4.4) 
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and using transformations ( [4.lD -( ^^ , Eqs. ( |3.5a| )-( ]3.5bD can be brought into forms analo- 
gous to the Dirac equation, where each particle appears as placed in the external potential 
created by the other particle. The wave equation satisfied by particle 1 becomes pT| : 



{ml 



m. 



2Pr 



2 2M 



+e 



-2Ua 



2x 



T2 



T T 
) .p 



-2ihe 



^ ™^ ' + 7il72lV2 + 715725 Vs + iI-iIUa + 



T T T T 

7/ .X 72 -x 



X 



.T2 



where we have defined: 
F I 



OF 



7aAt5 '^au 



(a = 1,2) . 



0, 
(4.5) 

(4.6) 



The wave equation satisfied by particle 2 can be obtained from Eq. (^4.5|) by the re- 
placements : pi ^ —P2, X ^ X, mi ^ m2, 7i ^ 72 • We recognize that the scalar 
potential Vi acts as a modification of the total mass M of the fermions through the 
change M — ^ Me^^i while {mf — m^) is kept fixed. The timelike vector potential V2 acts 
as a modification of the cm. total energy through the change Pie while 

[pIi^ — P2l) — i^i ~ ^2) is kept fixed. The spacelike potential f/4 changes the orbital 
angular momentum operator from L to Le ^ (in the cm. frame) and the combination 
U4 + T4 of the spacelike potentials changes the radial momentum operator from pr to 
(in the classical limit). The pseudoscalar potential appears only in spin- 
and ^-dependent terms. 

In QED, the summation of leading infra-red effects of multiphoton exchange diagrams, 



as described above |^, leads to the following expressions for the timelike (V2) and space- 
like (t/4 and T4) parts of the electromagnetic potential in the Feynman gauge: 



V2 



- In iH 

4 V Plt 



Vo 



. 



(4.7) 
(4.8) 



Potentials ( [4.7|) -( ^78D are compatible with the minimal substitution rules proposed long 
ago by Todorov for spin-0 particles, on the basis of an identification of the two-particle 
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motion in the cm. frame to that of a fictitious particle with appropriately defined reduced 



mass and energy |TH]. These rules were extended to the fermionic case by Crater and Van 



Alstine |jT^. The above potentials were shown to reproduce the correct O(a^) effects in 



muonium and positronium spectra |2^, |2T 



Similar results as above can also be derived with scalar photons contributing to the 
scalar potential Vi [Eq. ||TB|, ^ . 



The application of the previous summation method of the Feynman diagrams to the 
intercations of bosons shows, as one naturally expects, that the classical parts of the 
potentials in the fermionic case (written in the Pauli-Schrodinger form |^) and in the 



bosonic case coincide [^. Therefore, one can use unified potentials for both cases. 

It is not straightforward to generalize the above evaluation and summation techniques 
to other covariant gauges than the Feynman gauge. The presence of the additional gauge 
piece in the photon propagator breaks the permutational symmetry used in the previous 
calculations and renders their evaluation rather complicated. This is why these potentials 
will be evaluated from the Feynman gauge, using the infinitesimal transformation laws 
obtained in Sec. 3. 

To have a rough idea of the expected results, we first consider in some detail the 
one-photon exchange approximation . The corresponding potential is then local, with- 



out further approximation. Indeed, Eq. (|3.10|) , specialized to the one-photon exchange 



diagram, projects, with the constraint condition (|3.2| ), the photon propagator in mo- 
mentum space on the surface = 0; the potential in x-space is obtained with the 
three-dimensional Fourier transformation with respect to k'^. One finds for the photon 
propagator in three-dimensional x-space, in the gauge ^, the expression: 

D,Ax^) = iz { C + Cii - 1) + 1 ; (4.9) 



it yields the following expressions for the potentials in the gauge ^: 

[The dot operation is defined in Eq. ([4.6|) .] To this order, is independent of ^. 

After replacing these potentials in the wave equation ([4.5| ) (and in the equivalent one 
of particle 2) and designating by \E'^ the corresponding wave function, it can be seen that 
the wave function transformation 
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= -{Fx'.p' +p\x'F) 



(4.11) 



2' " " ' ' Air 2PLr ' 

removes, to first order in a, all the ^-dependent terms from tlie wave equation and gives 
back tlie wave equation in tlie Feynman gauge. Tlie operator Sq can also be written in 
tlie following form: 



Sn 



I 



Hn 



2Pr 



Hn,A(x^ 



(4.12) 



where A{x^) is the three-dimensionally reduced expression of the gauge propagator A(x) 
[Eq. ( p.3| )] (calculated by the Schwinger parametrization and dimensional regularization): 

Jk'^.x'^ -[ 

(4.13) 



A(x' 



1 



(27r)3 (k^^ + ie)^ 

The above study can also be repeated for the bosonic case, the same results as in Eqs. 
([4.11|) - (|4.12|) being found. [The wave equations for vector interactions with bosons can 
be found in Refs. |16|, [1^. In the first paper of Ref. |T^, the eigenvalue equation in the 
Feynman gauge is given by Eq. (5.12) with the identifications {1 — B) = {1 — A)^^ = 

To investigate the transformation laws for the higher order diagrams, we go back to 
the general case of Eqs. ( p. 131) and ( p.l7| ). We use the wave equation (|3^ ), together with 
relation ( |3.1CI| ), in its integral form, valid for a bound state, 

^ + gof{l-gof)-'^ = , (4.14) 

and add it, multiplied with an appropriate (nonsingular) factor, to Eq. ( p.l3| ). We obtain: 



6^^ 



-le 



'^iH^ {(AGoT) + (A^oT)} (1 - g^T)-^^ - -l-^e^^^i^oA^ . (4.15) 

The quantity (A^o^") is the constraint diagram counterpart of the amplitude (ACqT) ||20| 
and the integration inside it is three-dimensional, after constraint ( p.2|) is used. 

We can still add to Eq. (|4.15| ) the contribution of the crossed diagram counterpart of 
(AGqT), denoted by (AGoT)cr, in which the gauge propagator A crosses the scattering 
amplitude T (see Fig. 2). This is possible since (AGo^)cr does not have a pole at the 
bound state position in the s-channel and hence one can apply the operator (1 — g^T)'^ 
on \l/ and obtain zero. Thus, Eq. (|4.15| ) becomes: 



2P, 



-le 



'6^Ho {(AGoT) + (AGoT)„, + (AgoT)} (1 - goT)-'^ 



2P, 



-ie^^iHoAm . 



(4.16) 
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The sum of the amphtudes {AGqT), {AGoT)cr and (AgoT) can be evaluated at leading 
order of the infra-red counting rules of QED with the eikonal approximation p5|, |26|, |27 



adapted to the bound state problem. This approximation was verified to yield the correct 
results for the leading terms of the two-photon exchange diagrams and then was general- 
ized to higher-order diagrams [^. It consists of making the following approximations in 
the fermion propagators: 



Gi{pi - ki) 

G2{-{p2 + k2)) 



—2pi.ki + ie 

i 



illLPlL + mi) - 7iLfclL 
-12LP2L +1^2) - 72L^"2L 



(4.17) 



2^2-^2 + ie 

and of neglecting, at intermediate stages of the calculation, momentum transfers relative 
to subgraphs of a given graph. Neglecting thus in the amplitude T, in Eq. ( [4.1(j| ), 
the momentum transfer, the calculation becomes similar to that of two-photon exchange 
diagrams. (For positivity reasons, we also retain the quadratic term k1 in the product 
GiG2.) One finds: 

Ho 



(AGoT) + (AGoT),, + (A^oT) - A(2 



)T 



(4.18) 



4:PilP2L ' 

At leading order, the term Hq/ {^PilP2l) is equivalent to —^0"^; neglecting quantum effects, 
the latter can be brought to the utmost right and replaced there by V . One finally obtains: 



5^^ 



2Pi 



(4.19) 



where Viead,i is the infra-red leading part of V (in the gauge ^), i.e., the timelike component 
of the vector potential and where 7il72L is replaced by —1: 



Viead,i = -tanhV25 . 

(In the Feynman gauge V2 is given by Eq. (^4.7|) .) 
Similarly, Eq. (|3.171) yields: 



(4.20) 



M9o' + V)]^ 



2P. 



-te 



(4.21) 



To integrate Eq. ( [4.19| ) up to finite ^'s, we bring the operator Hq to the utmost right 
and use the equation of motion (with the approximation Hq ^ —^PilP2L%^)'- 



= 



2Pr 



ie^5i\^ Ho, A(l + Viead,if + ApilP2lA{1 + V^w,?) Vw,?} ^ • (4.22) 



The solution of this equation is: 



^6 = f^(6,6)% , 



2Pr 



Ho, A(l + Vlead,^)^ + 4pilP2lA{1 + Vlead,^fVlead,^'j , (4.23) 



where V is the path ordering operator. 

In the following, we shall study an approximate form of this transformation law. To 
this end, we adopt two simplifications. First, we assume that in the commutator, in W{C,), 
the potential Viead,^ can be approximated by its expression of the Feynman gauge: 



Vlead,^ ~ Vi,ad,F = \ ' • (4.24) 

Second, we assume that the two operators in W{^) are commuting objects. With these 
approximations the gauge transformation operator U{^2,^i) takes the form: 



f/(6,ei) ^ T(e2,ei)f/(6-ei) , 



UiO = e'^ , S 



2Pl 

2 



Ho,A{l + Viead,F) 



T(6, = exp 4piz.P2lA(1 + l^ead,c)'^ead,c} • (4.25) 

(T(^25'Ci) and f/(0 are supposed to be commuting.) As we shall see in Sec. 5, the above 
approximate forms provide the main qualitative properties of the gauge transformations 
of wave functions and potentials. 

The wave equation operator {qq^ + V) transforms as: 

[^-1 + v]^ = T-i(e2,ei)f/(6 - ei)(^o"' + v)u^{^2 - 6)^-^(6,^0 • (4.26) 



Transformations ( [4.25| )-( p:.26| ) ensure the (weak) invariance of the norm ( |3.19| ). 



Among the two transformation operators U and T, it is the former which is the 
nontrivial one, generating local transformations in x-space, while the latter acts as a 
multiplicative factor. In the following we shall focus our attention on the properties of 
the operator U. 



19 



5 Properties of gauge transformations in the local 
approximation of potentials 



This section is devoted to the study of the properties of gauge transformations in the 
local approximation of potentials, implemented by the operator U [Eq. ( [4.25| )]. There 



is a complete similarity between the cases of bosons and fermions (the operator U is 
spin independent) and for this reason we shall concentrate on the case of fermions only. 
The generator S of the transformations will be written in a form similar to that of Eq. 
(|4.11|) , which is more tractable for practical calculations. Thus, the gauge transformation 
operator U{^) is defined as: 



/ = ^ ^[r(i + w)1 , W = TTTT^ I • (5-1) 



The operator U{S,) acts through changes of the variables and . We denote by 
x-^(^), r{S,) and p^{i) the new expressions obtained from , r and p^, respectively, after 
f/(^) has acted on them. We have: 



^liO = U{0^lU\0 , KO = U{OrU\0 , r = sj-x^^ , (5.2) 
from which we deduce the differential equations: 

^ = -mO>liO , ^ = -/(r(O)r(O . (5.3) 

We notice that the variable x^/r remains unchanged under the action of f/(^), which acts 
as a local dilatation operator in a;-space, and hence it is sufficient to study the variation 
of r. 

The solution of Eq. (|5.3|) is: 

ri^) dz 



< , (5.4) 



from which we also deduce: 



(5.5) 



dr ff{f) 

The action of U{$,) on the momentum operator is more involved. The operator 
p^{^) is no longer parallel to p^ and has components along x^. For reasons that will 
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become evident below, we parametrize p^{C,) by means of two functions Ug^ = Ug^{r{C,), r) 
and Tg^ = Tgg(r(^),r) as follows: 



pliO = u{i)plu\i) 



.-2Ug^jr 



-2U, 



9i 



-2T, 



X 



X 



T2 



X .p 



+ihx„ 



X 



T2 



. (5.6) 



[The dot operation is defined in Eq. ( ^.61 ).] The last term, proportional to ihx'^, is fixed by 
the hermiticity condition. From the definition of ^6 obtain the differential equation: 

dpliO 



+5thxliOfiriO) + 2ihx^\i)xl{i)f{T{i)) . 



(5.7) 



[The dot derivations are with respect to X"^^(^).] Use in both sides of Eq. (5/7) of 
parametrization ( p.6| ) leads to differential equations concerning the functions Ug^ and Tg^: 



dU, 



9i 



/(r(0) 



9? 



2x''{0f{r{0) . 



(5.8) 



[The terms proportional to ihx^ do not lead to new conditions.] Taking into account Eq. 
(|0| ) and the boundary condition p^{C, = 0) = p^ , the solutions of Eqs. (|5.8|) are: 



U, 



9i 



2 \ r 



2 I /(r) , 



(5.9) 



These also imply the relation: 



4r' 



.dU, 



9? 



Q^2 



e^^9i - 1 . 



(5.10) 



In the nonrelativistic limit one has the behaviors 



U, 



9i 



uT + o 



2Af"^« ^""\M^) ' 2M 
In this limit, Eq. ( ^.10| ) reduces to the relation: 

J.NR ^ ^^2<^^9i 



VM2 



(5.11) 



-9i 



Qj.2 



(5.12) 



We next study the action of the operator U on the wave equation operator [qq + V). 
We first consider the operator g^^ [Eq. (|3.6|) ], which is composed of the Dirac operators 
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(7i-Pi-F"^i) and (— 72-P2=F^2)- The Dirac operator (71.p1— mi), say, becomes (71.p1— mi)^, 
where only the operator has changed, according to the transformation law ( |5.6| ). The 
operator (71 .pi — mi)^ has the same structure as the wave equation operator ( |4.5| ), in 
which Vi = V2 = V3 = and f/4 = f/^^, T4 = Tg^. [The terms proportional to the matrices 
present in Eq. (|4.5|) , mutually cancel out when expressions (|5^ ) are used for Ug^ 
and Tg^.] A similar conclusion is also obtained with the Dirac operator (—72.^2 — ^2)£,- 
The two operators (71.^1)5 and (—72.^2)5 (strongly) commute and hence (^o'"'^)^ is a well 
defined operator, in which the ordering of the Dirac operators is irrelevant. 

The action of the operator U on the potential V is obtained by the replacement in it 
of r by r(^), according to the transformations (|5.2| ) and (|5.4|) : 

% = V{r{i)) . (5.13) 

Examining then the norm of the wave function [Eq. (|3.19| )] (in the kernel of which, 
after the evaluation of the action of one uses the equations of motion), one deduces 
that the passage to the wave function characterized by a norm of the type ( [4.3|) , is 
again obtained with transformations of the type (|4.1| )- (|4.2| ): 

% = tanhl^^, ^5 = (coshF^)^^. (5.14) 

The Dirac type wave equations satisfied by have the same structure as Eqs. ( [4.5|) , in 
which, however, the potentials Ug^ and Tg^ have been added up to the existing potentials 
Ui{r{C,)) and T4(r(^)) of V^. This feature indicates us that these wave equations could 
also have been obtained from the following wave equation satisfied by a wave function 
defined below: 

(^0"' + VI) % = 0, (5.15) 

where is defined as: 

VI = t&nhVl , (5.16) 
and has the following timelike (V^'g) and spacelike (f/4g, T^^) components: 

^4 = = r2(r(0) , 

U',^ = U,^ + Ug^ = UMO) + UgMO,r) , 

= T,^ + Tg^ = TMO)+TAr{0,r) . (5.17) 
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The wave function is related to by the transformation: 

% = (coshV^/)^^ . (5.18) 
The relationship between and is: 

^ = (cosh l^^') (cosh V^g) -1^5 . (5.19) 
Therefore, the two wave equation operators 1%^ + V]^ and {gQ^ + V^) are equivalent: 

[%' + ^ + VI ■ (5-20) 

The advantage of the representation \E'^ is that its wave equation operator has the con- 
ventional form ( p.7| ) with a potential which is local. In this representation, among the 
three potentials V2, f/4 and T4, only V2 has a form invariant transformation law. If the 
scalar and pseudoscalar potentials, Vi and V3, were present, they would transform as V2. 
(In this case Viead [Eqs. (|5.1|) ] should also contain the scalar potential Vi.) 

The above transformation laws satisfy the group property, as can be checked by com- 
posing two succesive transformations; hence, they can be used starting from any gauge. 
Let us now return to the explicit expression of the function /, Eq. (p.l|). Equation 
^) then yields: 



1 / 2 \ 1 / 2 A L 2a 

exp — — = exp t, ] , X 



(2x(0 - 1) \x{0 - 1/ (2x - 1) Vx - 1 ^/ ' V Plt ' 

(5.21) 

It does not seem possible to express r(^) in compact form in terms of r and ^. However, 
the above equation provides easily the asymptotic behaviors of r(^): 

^iO = r-^ + o(^) , (5.22a) 



r->oo 2Pl \r , 



(5.22b) 



We have plotted, in Fig. 3, the curves r(^) (in units of 2a/ Pl) for three values of the 
gauge parameter, ^ = — 2 (Yennie gauge), ^ = (Feynman gauge) and ^ = 1 (Landau 
gauge). It is observed that the large-distance behavior ( ^.22a| ) is reached very rapidly. 
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From Eqs. ( |5.1| ) or Eq. ( |5.5| ) one also obtains the asymptotic behaviors of f{r{C,))/ f{r): 



f{r) ' 



1 + 



+ 0{- 



(5.23a) 



1 , 1 
/(r) r^o 2 V 2a 



1 + 



+ 4(l-e~^l 



2q: 



Oir) . 



(5.23b) 



Equations ( |5.22a| )-( ^.23b|) , together with Eqs. (|5.9|) and ( p.l7|) , yield the behaviors of 
the new potentials in the corresponding limits. The large-distance expansions ( |5.22a|) and 
(|5.23a|) are particularly relevant in the nonrelativistic limit. 

Of particular interest is the short- distance behavior of the effective potential. Ac- 
cording to Eqs. (|5.1| )- (|5^ ), the wave function \l/^(x"^) is equal to '^{x^{^)); Eq. ( |5.22b|) 
indicates us that r(^) behaves like r near the origin and therefore the behavior of the 
wave function does not change there under the gauge transformation; this in turn means 
that the dominant short-distance singularity of the effective potential has remained the 
same. These features can also be verified explicitly from the wave equations. 

To this end, let us consider, for the electromagnetic interaction, the Pauli-Schrodinger 
type eigenvalue equation obtained from the wave equation ([4.5|) (in the cm. frame) pT| : 



.4(f/4 + T4) 



2 _ 



rrin 



4P2 



^^^^ 



1) + 







(5.24) 



where 03 is a reduced (four- component) wave function and L is the orbital angular mo- 
mentum operator; the dots stand for spin- and /i-dependent terms, which are not relevant 
for the present purpose. The dominant short-distance singularity is provided by the term 
g4(V2 + f/4 + T4)_ ^YiQ Feynman gauge, [Eqs. ( [4.7|) - (f4.8|) ], it yields the singularity 
a'^/r'^, which is attractive and produces at the critical value = | the fall to the center 
phenomenon. A detailed analysis of this problem has shown that the theory undergoes 
at this value of a a chiral phase transition |2^; such a conclusion has also been reached 
from the Bethe-Salpeter equation in the ladder approximation [^] and from lattice theory 
calculations [^, Considering Eq. ( |5.24|) in the gauge ^, the potentials V2, U4 and 
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Ti become replaced by ^'5, and T'^^^ respectively. We find that the modification of 
the coefficient of the short-distance singularity coming from the form invariant part of 

+ ^4g! from + is cancelled by that coming from the form noninvariant 
part Ug^ + Tgg, and therefore the same singularity cP' jr^ emerges again. 

One consequence of this result is that the critical coupling constant has the same 
value \ in all gauges. This is a consistency check of the formalism, since a, representing 
here the invariant charge, should lead to a gauge invariant critical value etc. (In the present 
approximation, where radiative corrections are neglected, the physical and bare charges 
are identical; in any event, both quantities should be gauge invariant.) This is in contrast 
with the results obtained from the Bethe-Salpeter equation in the ladder approximation, 
where the value of is gauge dependent (equal to 7r/4 in the Feynman gauge and to vr/S 



in the Landau gauge |]29|, |32|). 

We also emphasize the particular role played by Viead^ in / [Eq- ( 101) ], in the short- 
distance behavior of the potentials in the gauge ^. If V^ead were absent, then the exact 
solution of Eq. ( p.4| ) would be r(^) = r — a^ji^Pj)) (the same as the asymptotic behavior 
(|5.22a|) ), producing in the Coulomb potential a singularity shifted to the position r = 
a^/{2Pi). On the other hand, when Viead is present in /, the function rf (cf. Eq. (|5.4|)) 
vanishes like r when r tends to and, as a result, the singularity of the Coulomb potential 
becomes maintained at the position r = in the gauge ^. 

The results obtained so far with the approximations ( [4.25[ ) are not qualitatively mod- 
ified when the dependence on ^ of Viead,i is introduced with an iterative treatment. This 
is a consequence of the asymptotic behaviors (|5.22a|) - (|5.22b| ) and of the fact that the 
corresponding behaviors of the function A(l -|- VJead,^)^ in / [Eq. (|0|) ] are not changed. 
Similarly, one can also estimate the effect of the operator T [Eq. ( [4.25| )]. The integrand in 
the argument of the exponential in T is a positive function. The behavior of the integral 
can be studied in the asymptotic regions with the aid of relations (|5l22a|) - (|5.22b|) . In the 
limit r — 0, the argument of the exponential tends to zero, while for r 00, it tends 
to the constant value (^2 — ^i)q^^PilP2l/(2-P|) and generally remains a smooth function 
between these two limits. These properties justify the factorization approximation made 
at the level of the gauge transformation operator U [Eqs. ( |4.23| ) and ([4.25|) ]. 

Let us ffnally remark that transformations ( ^.1|) , because of their spin-independent 
character, can also be applied in the case of bosonic wave equations, where now the 
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operator Qq is equal to —Hq, and the distinction between the representations and 
becomes irrelevant. 
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6 Summary and concluding remarks 



Using connection with quantum field theory, we estabhshed the infinitesimal covariant 
abelian gauge transformation laws of constraint theory two-particle wave functions and 
potentials and showed weak invariance of the coresponding wave equations. Contrary 
to the four- dimensional case and because of the three-dimensional projection operation, 
these transformation laws are interaction dependent. 

Simplifications occur when one sticks to local potentials, which result, in each formal 
order of perturbation theory, from the infra-red leading effects of multiphoton exchange 
diagrams. In this case, the finite gauge transformation can explicitly be represented, with 
a suitable approximation and up to a multiplicative factor, by a momentum dependent 
unitary operator that acts in x-space as a local dilatation operator. The latter acts on the 
potentials through two kinds of modification: a change of the argument r of the potentials 
into a function r(^) and a functional change of certain parts of the spacelike components 
of the electromagnetic potential. The function r(^) is essentially dominated by its large- 
distance behavior, in which r is simply shifted by a constant value. The knowledge of these 
modifications allows one to reconstruct, starting from the Feynman gauge, the potentials 
in other covariant gauges. It was shown that the dominant short-distance singularity of 
the effective potential of the Pauli-Schrodinger type eigenvalue equation is gauge invariant 
with a critical vahie of the coupling constant equal to 1/2. 

The above results allow the search for optimal gauges when incorporating new contri- 
butions into the potentials, which might come either from QED, or from other interactions. 
For instance, it is known that vacuum polarization affects only the transverse part of the 
photon propagataor. Therefore, the choice of the Landau gauge for the introduction of 
the effective charge seems to be most indicated. 

The similarity in structure between QED and perturbative QCD, up to the color 
gauge group matrices and the difference in the effective charges, allows us to envisage the 
consideration of many of the previous results in problems of quarkonium spectroscopy, 
where one has also to incorporate at large distances the effects coming from the confining 
potential. Here also the search for optimal gauges may become useful for subsequent 
applications. 
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Figures 



Fig. 1. Infinitesimal change, under gauge transformations, of two-particle Green's func- 
tion. (One-particle radiative corrections are neglected.) The dashed line, including its 

vertices, is the gauge propagator ie^d^/ {k'^Y . 

Fig. 2. The three diagrams entering in the evaluation of the interaction dependent part 
of the change of the constraint theory wave function. The shaded box is the off-mass shell 
scattering amplitude. The cross indicates the constraint diagram. 

Fig. 3. The curves r(^) versus r (in units of lajPi) for three values of the gauge 
parameter, ^ = — 2 (Yennie gauge), ^ = (Feynman gauge) and ^ = 1 (Landau gauge). 
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